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Most approximate exchange-correlation functionals used within density functional theory are con- 
structed as the sum of two distinct contributions for exchange and correlation. Separating the 
exchange component from the entire functional is useful since, for exchange, exact relations exist 
under uniform density scaling and spin scaling. In the past, accurate exchange-correlation potentials 
have been generated from essentially exact densities but they have not been correctly decomposed 
into their separate exchange and correlation components (except for two-electron systems). Using 
a recently proposed method, equivalent to the solution of an optimized effective potential problem 
with the corresponding orbitals replaced by the exact Kohn-Sham orbitals, we obtain the separation 
according to the density functional theory definition. We compare the results for the Ne and Be 
^ . atoms with those obtained by the previously used approximate separation scheme. 
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' Within density functional theory (DFT), the ground state energy of an interacting system of electrons in an external 
potential can be written as a functional of the ground state electronic density In the Kohn-Sham formulation of 
density functional theory [Q] , the ground state density is obtained as the density of a system of non-interacting electrons 
in an effective local potential. Although density functional theory is in principle exact, the energy functional contains 
an unknown quantity, called the exchange-correlation energy, E xc [p]. The effective potential for the fictitious non- 
interacting system is the sum of the external potential, the Hartree potential and the exchange-correlation potential, 
which is the functional derivative with respect to the density of E xc [p]. The density functional theory definition 
of the separate exchange and correlation components of E xc [p] is based on the non-interacting system and is such 
that the resulting exchange functional has properties that are useful guides in the construction of an approximate 
exchange. Consequently, most approximate exchange-correlation functionals are also constructed as the sum of two 
distinct contributions for exchange and correlation. 

In the past, exchange-correlation potentials and energies, of varying degrees of accuracy, have been determined by 
generating a density for the system of interest and then computing an exchange-correlation potential that yields the 
desired density as the ground state solution for the fictitious non-interacting system. In this context, researchers 
have used charge densities calculated by quantum chemistry methods for atoms [p|-[l3[ and molecules |14|dTq|, as well 
as Quantum Monte Carlo methods for atoms jl7|,[lj| and for a model semiconductor [Q. The subsequent inverse 
problem, namely the search of the corresponding exchange-correlation potential, has been performed using a variety 
of different techniques. For example, the exchange-correlation potential has been determined by expanding it in a set 
of basis functions and varying the expansion coefficients to reproduce the accurate density ||,|l9| . 
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With the exception of two-electron systems, these accurate exchange-correlation potentials have never been separated 
into their exchange and correlation components according to the density functional theory definition. For many- 
electron systems, an approximate scheme was used where the exchange potential was defined as the potential yielding 
the Hartree-Fock density and the correlation potential as the difference of the accurate exchange-correlation potential 
and this approximate exchange potential. In this paper, following the approach proposed by Gorling and Levy p2[ , 
we obtain for the first time the correct separation of accurate exchange-correlation potentials for the Be atom and 
the Ne atom. 

In Sec. 0, we briefly introduce density functional theory and its Kohn-Sham formulation. In Sec. Ill, we derive the 
formulae used to determine the decomposition of the exchange-correlation potential into exchange and correlation. A 
comparison with approximate separation schemes is given in Sec. jtVl In Appendix [A], we describe the method for the 
special case of closed shell systems. 



II. THEORETICAL BACKGROUND 



Density functional theory provides an expression of the ground state energy of a system of interacting electrons in 
an external potential as a functional of the ground state electronic density Q| . Let us assume for simplicity that the 
spin polarization of the system of interest is identically zero. In the Kohn-Sham formulation of density functional 
theory (2| , the ground state density is written in terms of single-particle orbitals obeying the equations in atomic units 
(h = e = m = 1): 

|-^V 2 + t; cxt (r) + J ^Ldr' + v xc ([p}- ) v)U i = e i t/j i , (1) 

where 

p(r) = J2\Mr)\ 2 - (2) 

i=i 

The electronic density is constructed by summing over the JV lowest energy orbitals where ./Vis the number of electrons. 
fext(r) is the external potential. The exchange-correlation potential i; xc ([p];r) is the functional derivative of the 
exchange-correlation energy £" xc [p] that enters in the expression for the total energy of the system: 

E = -\ J ^V 2 ^ dr + J p (r) v ext (r) dr + ~ J J P ^ M ^ dr dr' + E xc [p] . (3) 

The exchange-correlation functional is written as the sum of two separate contributions for exchange and correlation, 

E xc [p] = E x [p] + E c [p] . (4) 

The definition of the exchange energy is in terms of the non-interacting wave function <I>o, the Slater determinant 
constructed from the Kohn-Sham orbitals, as 

E x [p] = <$ | V cc \*o)-\JJ ^^dr dr', (5) 

where V ee is the electron-electron interaction. This definition differs from the conventional quantum chemistry defi- 
nition of E x as the exchange energy in a Hartree-Fock calculation, given by the same expression as in Eq. but with 
the Kohn-Sham determinant replaced by the Hartree-Fock determinant. The separation of the exchange-correlation 
functional into exchange and correlation yields a corresponding splitting of the exchange-correlation potential into 
v x ([p] ; r) and v c ([p] ; r). In this formulation, the essential unknown quantity is the exchange-correlation energy E xc [p\. 
If the functional form of E xc [/?], and consequently the exchange-correlation potential, were available, we could solve 
the JV-electron problem by finding the solution of a set of single-particle equations. 

The exchange functional, as defined in Eq. [|, scales under uniform density scaling, p\(r) = A 3 p(Ar), Q as 
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E x [p A ] - XE X [p] , 



(6) 



and its spin polarized version is simply given in terms of the unpolarized exchange functional pjj as 

E x [PhPl] = \{E x [2pi]+E x [2p l ]}. (7) 



Clearly, the separation of the exchange-correlation functional according to Eq. |j| is useful since, for exchange, only an 
approximation for the unpolarized functional needs to be sought and the behavior under uniform scaling determines 
how derivatives of the density combine with the density in an approximate exchange functional: 

£ x approx [P\ = / P(r) 4/3 F{ |Vp(r)| Mr) 4 / 3 , V 2 p(v)/p(vf'\ . . . ) dr. (8) 



III. SEPARATION OF Vxc INTO Vx PLUS V c 



For the special case of two electrons in a singlet state, the separation of the exchange-correlation potential into 
exchange plus correlation is quite simple since the exchange potential is simply given by the condition that it cancels 
the self-interaction term in the Hartrcc potential. On the other hand, for many-electron systems, this decomposition 
into exchange and correlation components has never been done. In previous work ]5|,|6|,p|, pT|Jl^ ] , the exchange potential 
was defined as the difference of the effective Kohn-Sham potential yielding the Hartree-Fock density and the sum of 
the Hartree and the external potentials. The correlation potential was then obtained as the difference of the exchange- 
correlation potential corresponding to the exact density and the above approximate exchange potential. Note that this 
"exchange" potential is not an exchange-correlation potential since we are subtracting the wrong external potential: 
the Hartree-Fock density is the true ground state density for a Hamiltonian with an external potential different than 
the original one. However, it is also not the exchange potential corresponding to the Hartree-Fock density (although 
very close to it) since it is not the functional derivative with respect to the density of the exchange energy evaluated 
for the orbitals obtained from the effective potential yielding the Hartree-Fock density. Therefore, this separation 
scheme is incorrect: it involves two densities, the exact and the Hartre-Fock densities, and, moreover, the potential 
used for exchange is only approximately equal to the exchange potential corresponding to the Hartree-Fock density. 

We follow Gorling and Levy p2| ] in showing how to separate the exchange-correlation potential into exchange and 
correlation. We consider a spin unpolarized system. If we assume that the density p is non-interacting v-representable, 
it can be expressed as in Eq. || in terms of single-particle orbitals {ipi} of the Kohn-Sham potential v s (r), 

v s (r) = Wext (r) + f ^Ur' + i, xc {[p] ; r) . (9) 
J r - r' 



The exchange energy is a functional of the density but can also be expressed in terms of the Kohn-Sham orbitals {ipi} 
(Eq. g) as 

E * w = - 5 EE^,^ / / | r l r , drdr ' (10) 

i=l j=l ' 11 



where the (^-function is over the spin quantum numbers of the i-th and j-th orbitals. We evaluate the functional 
derivative of the exchange energy functional with respect to the Kohn-Sham potential as 



On the other hand, since the exchange functional can be written as a function of the orbitals (Eq. [n]), we also have 

SE x [p] " f f 5E x [p] SMr') 5E*[p] S^(r') \ 
Sv s (r) V^i(r') M') WiP) M r ) ) ' 
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If we combine Eqs. |ll] and [T^, we obtain the integral equation 

|« x ([p];r')X;(r',r)dr' = Q(r), (13) 

where the kernel /C(r',r) and the right hand side Q (r) depend on the orbital {ipi} and their functional derivative 
with respect to the potential v s (r) . This integral equation is equivalent to the one solved in the optimized effective 
potential method (OEP) where the KS orbitals are replaced by the OEP orbitals |23|] . The functional derivatives of 
the orbitals Sipi(r)/5v s (r') can be expressed in terms of the Green's function Gi (r, r') as 

|^ = -G i (r,r')M*')> (14) 
Sv s (r') 

where Gi (r, r') satisfies the differential equation 

iv 2 + v s (r) - e^j Gi (r, r') = 8 (r - r') - ^(r)^*(r'). (15) 

By knowing the exchange-correlation potential, the KS orbitals and eigenvalues, we can compute the Green's functions 
{Gi} and, consequently, the kernel K, and the function Q. If we express the exchange potential as a linear combination 
of basis functions, Eq. [l^ can be rewritten as a non-homogeneous set of linear equations for the coefficients of the 
expansion of the potential in the basis set. 

Once the exchange potential is determined, the correlation potential is simply obtained as the difference: 

u c([p];r) = v xc ([p];r) -v x ([p];r). (16) 

In Appendix [a| the equations derived in this section are rewritten for the case of a closed shell system. 

IV. COMPARISON WITH APPROXIMATE SEPARATION SCHEMES 

For the Be atom and the Ne atom, we calculate the exchange potentials as explained in the previous section. The 
correlation potentials are determined as the difference of the accurate exchange-correlation potentials and the exchange 
components (Eq. |l6|). We already mentioned that, in the past, an approximate "exchange" potential was instead used, 
given by the effective potential yielding the Hartree-Fock density minus the Hartree and the external potentials: 

v x ([pM,r)=v s ([pM,r)- [dr'^^l-v ext (r), (17) 



where we introduced an explicit dependence of the Kohn-Sham potential v s on the density reproduced by v s . The 
potential v x is very close to the exchange potential corresponding to the Hartree-Fock density, w x ([phf]; r). To 
determine i> x ([jOhf]; r), we can use the same scheme explained in the previous section with the orbitals given by the 
Kohn-Sham orbitals corresponding to the effective potential yielding the Hartree-Fock density instead of the exact 
density. We denote by and vf the correlation potentials determined as the difference of the accurate v xc and 
■M[PHF];r) and w x ([/5 H f]; r) respectively: 

v^{r) =v xa ([p];r) ~ v x ([p HF ];r), vf(r) = v xc ([p];r) - w x ([phf]; r). (18) 

As discussed in Ref. Iljj, for two-electron systems, £ x ([/9hf]; r) = Vx([phf}',y) and, consequently, v^(r) = vf(r). 
Further, it was empirically found that the difference between v c (r) and i^' B (r) is small on the scale of v c (r). 

Here, we find that even for the many-electron atoms Be and Ne, the difference between i; x ([p];r), i; x ([phf]; r) and 
"5 x ([/9hf]; r) are almost not visible on the scale of v x ([p];r). As shown in Fig. [j], the difference between v^(r) and 
v®(r) is just barely visible even on the more expanded scale of v c (r), This agreement is expected since the HF and 
the OEP densities are very close to each other and, for the OEP density, the agreement would be perfect. On the 
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other hand, the difference between either v c (r) and v£(r) or v c (r) and v^(r) is visible. For both atoms, the exact and 
approximate potentials are clearly different, although the shapes are very similar. The similarity of the exact and the 
approximate potentials justifies the use of the approximate scheme used in earlier work. 
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FIG. 1. Comparison of the correlation potentials of Be and Ne. v c is the correct correlation pot ential from Eq. hq, v c (r) and 
vf?(r) are the approximate correlation potentials constructed from the HF density using Eq. [18l (r) and vf(r) are nearly 
indistinguishable on the scale of these plots. For Ne, there is some uncertainty in the potentials for r < 0.4 ao. 
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APPENDIX A: CLOSED-SHELL CASE 

In Sec. |l], we presented the theory of separation of the exchange-correlation potential into the exchange and correlation 
components for spin-unpolarized systems. In the present section, we restrict ourselves to the special case of a closed 
shell atom. For a closed shell atom, the self-consistent solutions of the Kohn-Sham equations (Eqs. |l| and |J) can be 
factorized as the product of radial and angular components: 

&(r) - Ri(r)Y hmi (v) = M2y ( f) (A1) 

r 

where li is the angular momentum quantum number. Using this expression for the orbitals, the density (Eq. |J) can 
be rewritten as a sum over the occupied shells: 

i=l 

where N s is the number of occupied shells and fa is the occupation number of the i-th shell, /, = 2(2^ + 1 ). 
Following the derivation by Slater [M, we rewrite the exchange energy (Eq. |l(J) as 

N s k+lj 

EAp\ = -X] V /(2/ i + l)(2/ J - + l) c k (k,0;l J: 0)G k (n u l l ;n J ,l J ), (A3) 

id k=\U-lj\ 
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where the coefficients c k incorporate the integrals over 9 and are tabulated in Ref. p3j and G k is given by 

G k (n i ,l i ;n j ,l j ) = f dn f dr 2 fa (n ) fa (r 2 ) fa (r x ) & (r 2 ) 4fi > ( A4 ) 



with r< = min{ri,r , 2} and r> = max{?-i, r2}. 

The functional derivative of the exchange energy with respect to the effective Kohn-Sham potential (Eqs. [ll] and |T^ ) 
can here be obtained taking into account that the density depends only on the radial components of the Kohn-Sham 
orbitals. Eq. [ll] is therefore equivalent to 

SE X [ P ] „ [°° AJ ^ ^ ^ f ^Sfa(r') 
Sv s (r) 

while Eq. reduces to 



SExlp] f°° j , SEUlp] 8fa{r') 



Sv s (r) J j^f 8fa(r') 8v a {r) 



(A6) 



Eqs. [A5| and |A6| can be combined to give the following integral equation for the exchange potential: 

dr' v x (r')K(r',r) = Q(r). (A7) 



The functional derivative of the exchange energy (Eq. [AJ) with respect to the radial orbital is given by 



8fa(r) 



121- + 1 f°° r k 

^EV^TT c k (k,0;lj,0)fa(r) dr* fa(r 2 )fa(r 2 )^. 



(A8) 



Gi(r,r%(r'), (A9) 



The functional derivative of the radial orbital, 8fa(r)/Sv s (r'), is expressed in terms of the Green's function Gj(r, r') 
as 

5v„(r') 

where Gi{r,r') satisfies the following differential equation: 

W + ^ ^ + < r ) - e ] G ^ r ') = <^ ~ r ') ~ Unfair')- (AlO) 

This equation can also be derived by starting from Eq. [l5| in polar coordinates and projecting out the radial component. 
It can be easily checked that Gi has the following expression: 

G«(ry)= E (All) 
e? — £i 

where the sum is over all the orbitals, except the «-th one, with angular momentum quantum number ij. In solving the 
differential equation for Gi, we set r ^ r', divide by fair 1 ) and determine Xoutjf) and Xm(f) as solutions of outward 



(r < r') and inward (r > r') integration. 0(r) is a homogenous solution of Eq. AlO and can be added to Xout(r) and 
Xi n (r) as a on t(j>(r) and CKi n </3(r) respectively. The difference a ou t ~»in is determined by imposing continuity on Gi and 
the sum a out + ct'm by requiring that 

„'« r « r >l£H (A12) 



G 



which follows from the normalization of <fii(r). Finally, we obtain 



4 % = Gi(rS) Mr) I dr" ^(r")G,(r", /), (A13) 



5v s (r') 
where Gi(r, r') is 

G t (r, r') = 9{r r>) |xin(r)^(r') 2 - \ [%;„(/) - Xout (r')] ^(r)^(r')} 



0(r' - , ) { X out(r)Mr') 2 + \ [Xin(r') - Xout(r')] *(r)*(rO \ ( A 14) 
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